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ABSTRACT 

Recently, Suzaku has produced temperature and entropy profiles, along with 
profiles of gas density, gas fraction, and mass, for multiple galaxy clusters out 
to ~r2oo (~ virial radius). In this paper, we compare these novel X-ray obser- 
vations with results from N-body + hydrodynamic adaptive mesh refinement 
cosmological simulations using the Enzo code. There is excellent agreement 
in the temperature, density, and entropy profiles between a sample of 27 mostly 
substructure-free massive clusters in the simulated volume and the observed clus- 
ters. This supports our previous contention that clusters have universal outer 
temperature profiles. Furthermore, it appears that the simplest adiabatic gas 
physics used in these Enzo simulations is adequate to model the outer regions of 
these clusters without other mechanisms (e.g., non-gravitational heating, cooling, 
magnetic fields, or cosmic rays). However, the outskirts of these clusters are not 
in hydrostatic equilibrium. There is significant bulk flow and turbulence in the 
outer intracluster medium created by accretion from filaments. Thus, the gas is 
not fully supported by thermal pressure. The implications for mass estimation 
from X-ray data are discussed. 
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Introduction 



Galaxy clusters are unique and potentially powerful cosmological probes of the Uni- 
verse. They are the largest gravitationally bound objects, having grown hierarchically 
within the large-scale cosmic web. As such, clusters are an important source of info rma- 
tion about the components of t he universe in which they formed and evolved (e.g., IVoit 



20051 ; iBorgani fc Kravtsovi 120091 ). Gravity drives structure formation within an expanding 
universe, with regions of density higher than the average becoming gravitationally bound and 
decoupling from the expansion. Clusters probe the high density tail of this cosmic density 
field. The number density of clusters is strongly dependent upon the specific c osmological 
model, especially when viewed as a function of redshift (e.g., iRosati et al.ll2002l ). The mass 
function of clusters and its evolution are also sensitive to as (a parameter that quantifies the 
RMS density fluctuations on comoving scales of 8 h -1 Mpc). The potential to use the clus- 
ter mass function to measure cosmological parameters is challenging observationally because 
mass is not a direct observational quantity. Rather, X-ray luminosity, X-ray temperature, 
thermal SZE (Sunyaev-Zeldovich Effect) integrated Y(oc J pTdl), or weak lensing shear are 
the observables which must then be translated into an estimate of cluster mass. 

Over the past decade, X-ray observations from ROSAT, Chandra, XMM-Newton, and, 
more recently Suzaku, h ave begun to produce significant samples of clusters out to z ~ 1.3 



Allen et al.ll2008l ). With these samples, it is becoming possible to dist inguish between 



{e.g. 

cosm ological models using the evolution of the cluster mass f unction (e.g., IVikhlinin et al. 



20091 ) and the cluster gas fraction (e.g., iRapetti et al.ll2008l ). The primary limiting factor 



in such applications remains accurate conversion between X-ray observables such as X-ray 
luminosity and temperature to cluster mass. This is complicated by the nonlinear baryonic 
processes at the cores of clusters, including radiative coolin g and non-gravitational heating 



from supernovae and active g alactic nuclei (AGNs; e.g., iRuszkowski &: Begelmanl 12002 



Heinz et al. 



e.g., see 



1998 



Burns 



Skillman et al. 



1998 



) and possibly nonthermal p rocesses involving cosmic rays 



Churazov et al. 


2008; 


Xu et al. 


in scaling relations ( 


Motl et al. 



20051 ; iNagail 120061 ) , a s well as errors created by assum i ng that cluster gas is in hydrostati c 



equilibrium (e. g. , see lLau et al.ll2009l ; iBurns et al.ll2008uRasia et al.ll2006t iNagai et aU20071 ). 



Because of the nonlinear nature of structure formation, t he regions outside the cores 



do not necessarily behave in this same complex manner (e.g., iRoncarelli et al.l 120061 ) . The 
peripheries of clusters have lower gas and galaxy densities, and long cooling times comparable 
to the Hubble time, thus potentially making the thermodynamics simpler than in the cores. 
At the same time, the outer regions of clusters are closer to the sources of accretion from 
filaments and, therefore, the gas dynamics may be more complex. 
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The intracluster X-ray gas in galaxy clusters is often assumed to be in hydrostatic 
equilibrium which relates the gravitational potential (<&) to the gas pressure (P) and the gas 
density (p) such that 

VP 

V$ = -— . 1 



p 

Applying Gauss's Law to the gravitational potential in the above, the cumulative mass is 



M(< r) 



AtiG 



VP 

P 



dA, 



(2) 



where G is Newton's gravitational constant and the integral is over a spherical surface area 
with radius r. If the cluster is further assumed to be spherically symmetric and the pressure 
arises only from thermal motions (P = pkT/[pm p ]), then 



M(< r) 



r 2 k 



pGpm 



p i 



dr dr 



(3) 



As seen in Equation (3), we need to measure both the gas densities and temperatures, along 
with their gradients, to calculate the hydrostatic equilibrium masses. Gas density profiles 
can be accurately determined from the X-ray surface brightness, Sx, since Sx oc p 2 with 
only a weak dependence on temperature. 

Temperature profiles, on the other hand, are more challenging since they require spatially- 
resolved X-ray spectroscopic observations acquired from relatively long integrations. This 
has been difficult, particularly to measure cluster temperatures beyond ~ 0.5r 2 o(jj Because 
of the low and stable particle background levels at its orbit, Suzaku has begun to change the 
landscape by producing X-ray temperature p rofiles out to ~r 2 nn for a small nu mber of clusters 



i ncluding PKS 074 5-171 (IGeorge et al 



20091) . Abell 1795 flBautz et al. 



200 91), Abell 399/40 1 



flFuiita et alJl2008f l. Abell 2204 flReiprich et al.ll2009h . and Abell 1413 flHoshino et al.ll2010h . 

In each cluster, the temperature is observed to decline by a factor of ~ 3 from the peak near 
the cluster core to regions at ~ T2oo- 

Using the Eulerian adaptive mesh refinement cosmology code Enzo for a ACDM uni- 
verse , our group first proposed a universal temperature profile for galaxy clusters (ILoken et al. 
20021 ). This average profile was well-fit by a power-law out to the virial radius. This pro- 
file agre ed well with the X-ray data available at that time for nearby clusters from Bep- 
poSAX (IDe Grandi fc Molendil 120021 ). Subsequent Chandra observations also appear con- 
sistent with such a universal temperature profile beyond the dense central cores for galaxy 



4 ?"2no is the radius at which the density is 200 times critical and 
|Brvan fc Normarj[l99l ) 



virial radius for a flat ACDM Universe 
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clusters (e.g., IVikhlinin et al.ll2006l ). However, these observed X-ray temperature profiles 
extended out only to ~ 0.5r 2 oo so it was unclear if the universal temperature profiles for 
numerical clusters are in agreement with the outer profiles for real clusters. 

In this paper, we compare the temperature, density, and entropy profiles for a new 
sample of numerical clusters generated using the Enzo cosmology code with new observations 
from Suzaku. We ask the question: Do real clusters follow a universal temperature profile 
out to the virial radius, as is predicted by numerical simulations with simple adiabatic gas 
physics? We also explore the implications of this particular form of the universal temperature 
profile for hydrostatic equilibrium, intracluster gas dynamics generated by on-going accretion 
in the outer periphery of clusters, and cluster mass estimation. 

In Section HI we describe the Enzo cosmology code and the numerical simulations used 
for the comparison with X-ray observations. In Section [31 we compare gas densities, temper- 
atures, and entropy profiles between our numerical simulations and observed clusters. Then, 
in Section HJ we explore the implications of the form of the observed and simulated cluster 
temperature profiles on hydrostatic equilibrium and cluster mass determinations. We end 
with a summary and conclusions in Section 



2. Simulations 



We use the Enzo adaptive me sh refinement (AMR) N-body + hydrodynamics cosmology 
code fervan fc Norman! Il997al lbl: iNorman fc BrvarJ Il999l ; b'Shea et aD |2004[ lO'Shea et al 
20051 ) to simulate a comoving volume of (128 h~ l Mpc) 3 with 256 3 root-grid cells and up 
to 5 levels of additional refinement. This simulation utilizes the ZEUS finite-difference 



method (jstone fc Norman 1992a, b). The AMR is controlled by refining any region that 
is overdense by a factor of 8 in either the dark matter or gas density. Th e peak resolution 
i s 15.6 /i _1 kpc (comoving). The initial conditions are generated from an lEisenstein fc Hu 
( 119991 ) power spectrum with a primordial spectral index n s = 0.97. We use the cosmolog- 
ical parameters (/i,0 M ,fi B ,0 CM/ ,fi A ,(T 8 ) = (0.704, 0.268, 0.0441, 0.2239, 0.732, 0.9), with 
h = iTo/(100 km s" 1 Mpc~ 1 ). The simulation was initialized at z = 99 and run until z = 0. 
The dark matter mass resolution is 3.12 x 10 9 /i~ 1 M . For a further expos i tion o f Enzo and 
its use in studying the statistical properties of clusters, see ISkillman et al.l ( 120081 ). 



For all the analyses of the simulations and observations throughout this paper, we will 
assume H = 70 km/sec/Mpc and the other cosmological parameters given above. 
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Table 1: Comparison of Properties of Simulated with Observed Clusters 



Cluster 


z 


M 200 (M ) 


r 200 (Mpc) 


PKS 0745-191 


0.103 


6.4 x 10 14 


1.72 


Abell 1795 


0.063 


8.6 x 10 14 


1.9 


Simulated Clusters 


0.05 


1.7- 10 x 10 14 


1.6-2.8 



3. Comparison of Simulated and Observed Clusters 



We selected the 40 most massive clusters in our numerical volume at z = 0.05 for 
comparison to the Suzaku observations because these clusters best match those observed 
(Tabled]). This sample was then edited by examining the 3D distribution of the dark matter. 
Any cluster with a secondary dark matter density peak located outside of 0.1r 2 oo was marked 
as "disturbed". Upon visual inspection, these clusters were confirmed to be ones with major 
substructure. This produced a "clean" sample of the 27 most massive clusters that are 
relatively substructure- free. However, a note of caution is warranted, as there is accretion of 
smaller halos (subclusters) present in nearly all clusters; such small halo accretion is better 
seen in temperature images where shocks from supersonically merging subclusters are more 
obvious (see Fig. [I]). 

A grid of images of gas densities, temperatures, and 0.5 — 12.0 keV X-ray emission is 
portrayed in Figure [TJ For each te mperature a n d den sity projection, we show a correspond- 
ing X-ray image. The method of ISmith et al.l ( 120081 ) was used for the radiative cooling to 
calculate frequency- dependent emission using the emission function from the Cloudy code 
(IFerland et al.lll998l ). Given a temperature and density of the gas, and assuming a metal- 
licity of 0.3 solar, this returns the total X-ray emissivity. The X-ray surface brightness is 
the integral along the line of sight of the X-ray emissivity. These images in Figure [1] show 
a good deal of non-circularly symmetric structures, often produced by accre tion along fil- 



aments, as also appear s to be seen in the clusters observed by Suzaku (e.g., iGeorge et al. 



2009: Bautz et al 



20091 ) 
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Fig. 1. — Grid of images of the 40 most massive clusters in the numerical sample selected 
from a cosmological volume at z = 0.05. The images are ordered by M 2 oo from the upper left 
(10 15 M ) to the lower right (1.7 x 10 14 M Q ) for each of the three panels of gas density (top), 
temperature (middle), and X-ray surface brightness (bottom). Field of view of each image 
is 8 h~ l Mpc. The range of properties for these clusters is listed in Table 1. This sample 
was reduced to 27 clusters for statistical analysis by eliminating obvious merger systems. 



3.1. Analysis Pipeline 



3.1.1. 3D Profiles 



For each cluster, we create spherically averaged radial profiles centered on the center 
of mass of the cluster. We track the following key fields: radius (in Mpc), overdensity of 
a sphere out to the current radius with respect to the background density (Q m p cr ), dark 
matter density, gas density, entropy, pressure, specific radial kinetic energy, specific kinetic 
energy, temperature, thermal energy, X-ray emission, total energy, and enclosed baryonic + 
dark matter mass. 

An important point in our profiling procedure involves the various options for weighting 
the profiles. In this study, we use two different methods. First we can use the simple method 
of weighting by the mass of the cell. Therefore, regions in radial shells in high-density regions 
will be weighted higher than those in underdense regions. Second, we weight by the total 
X-ray emission emanating from each cell. This biases towards overdense regions even more 
strongly since the X-ray emissivity scales roughly with p 2 . In Figures 3 and 4, we use X- 
ray weighting to compare with X-ray observables. In Figure 5, we use mass weighting to 
investigate the physical causes of deviations from hydrostatic equilbrium. 



Because observations yield inherently 2D images, it is important to compare projections 
of our clusters to the true observations. We do this by first creating halo projections of each 
cluster that encompass an 8 Mpc cube around the individual cluster. For each of these 
projections, we again have a choice between weighting by the density or X-ray emission. 
Unless otherwise noted, we will use only the projections weighted by the X-ray emission. 
From each of the 2D projections, we then seek to create radial profiles. Here we give details 
of creating 2D profiles of X-ray weighted temperature, and use the same technique in all 
other quantities. 

For each cluster, we first make a projection of the X-ray emissivity, ex, which yields the 
X-ray surface brightness, 



3.1.2. Profiles of 2D Projections 




(4) 



We then make a projection of the temperature, weighted by X-ray emissivity, 
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From these quantities, we create the profiles of the 2D projections by weighting the temper- 
ature projection by the surface brightness, 



where r and 9 are the normal polar coordinates. This results in a correctly weighted radial 
profile of X-ray weighted temperature. Note that this is different than simply creating a 2D 
profile from the projection of X-ray weighted temperature. We can compare this final step 
to the equivalent step in the creation of the spherically averaged radial profile of the same 
quantity using normal spherical coordinates, 




J Tp ro j : X-rayweighted(x, y)Sx,Proj(x, y)rd6 

J S x ,p ro j(x,y)rd9 



(6) 




/ j T(x,y, z)ex(x,y, z)r 2 sin9d9d(f> 



(7) 



f f €x(x,y, z)r 2 sin9d9d(fi 



For an illustration of this pipeline, see Figure [2j 
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Fig. 2. — Analysis Pipeline to construct both a 2D and 3D X-ray weighted profile of a halo. 
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3.2. Results 

Applying the approach described in Section [37T1 to the analysis of the simulated clusters, 
we produced radial profiles of density, temperature, and entropy as shown in Fig. [3j Here, 
we define entropy as K = kT(— £— )~ 2 / 3 . The profiles were constructed by binning in circular 
annuli centered on the peak X-ray emission. The temperature profiles were each normalized 
to their X-ray weighted average values between 0.2 to 1.0 r 20 o so we could easily compare 
clusters of different masses. Such normalization is less effective for density and entropy 
which were, instead, normalized by values at 0.1 r 2 oo arid 0.3 r 2 oo, respectively. The shaded 
regions illustrate the standard deviations in quantities derived from the sample of numerical 
clusters. This scatter generally grows with radial distance, especially for the density and 
entropy, beyond ~ 0.7r 20 o- This scatter reflects real variations from cluster to cluster in the 
outer regions produced by non-spherically symmetric accretion along filaments as mentioned 
above. The flat outer profile (> 0.5r2oo) seen for the average entropy of the numerical cluster 
sample in Figure [3] is, in part, due to the averaging of individual profiles, some with positive 
and some with negative slopes. 

The average temperature profile shown in Figure [3] is well- fit by a function of the form 



rp / \ n -3.2±0.4 



T 

-L avg 



;i.74±0.03) 



(0.64 ±0.10) 



r 200 



(8) 



As we discuss in Section HI the particular shape/slope of this power- law for the temperature 
profile has implications for understanding the gas dynamics and ICM pressure support in 
the outer regions of clusters. 

In Fig. [31 we also compare the numerical profiles with observed ones for two nearby 
clusters (PKS 0745-191 and A1795) with recent Suzaku X-ray measurements that extend out 
to ~ r 2 oo- A summary of the characteristics of these two observed clusters, along with the 
average properties of the numerical clusters, is given in Table [TJ Both observed clusters have 
small cool cores. Since our simulations were purposely constructed with simple adiabatic 
gas physics, we did not expect to match the central regions of these clusters. Instead, we 
focus on the outer profiles, beyond > 0.5r 200 . Other clusters recently observed by Suzaku 
have similar outer temp erature profiles with a factor of & 3 decli ne in temperature between 



the core and r 2 oo (e.g., iReiprich et al.ll2009t iHoshino et al.ll2010l ). 





Fig. 3. — 3D and 2D X-ray weighted average profiles, as described in Section 13.11 of nor- 
malized density, temperature, and entropy for simulated clusters. The shaded regions are 
the standard deviations of each quantity produced from the scatter within the sample of 
numerical clusters. The Suzaku data points for PKS 0745-191 and A1795 (measured for two 
separate analysis regions, north (N) and south (S) of the core) are shown for comparison. 
The blue and green shaded regions in the bottom entropy profile are for the N and S analysis 
regions of A1795. The temperature is normalized by the average value between 0.2 to 1.0 
r 20 o, whereas the density and entropy are normalized by values at 0.1 r 2 oo and 0.3 r 20 o, re- 
spectively. This difference in normalization produces a slight shift in the normalized entropy 
between simulations and PKS 0745-191. 
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It is clear from this figure that the outer cluster radial distributions of density, tempera- 
ture, and entropy for PKS 0745-191 and A1795 agree well with the simulated clusters. From 
Figure [31 we draw three conclusions about the regions outside the cluster cores. First, the 
observed clusters follow "universal" temperature, density, and entropy profiles as character- 
ized by numerical simulations out to the outer bounds of the clusters (r 20 o). Second, the 
simplest intracluster medium (ICM) gas physics with an adiabatic equation of state is suffi- 
cient to characterize the thermodynamic gas properties of the peripheries of these clusters. 
With the present Suzaku X-ray data, added gas physics such as cooling, non-gravitational 
heating ( via e.g., low powe r radio or X-ray AGNs associated with non-central cluster galax- 



ies [e.g. 



Hart et al.l 120091)) or nonthermal pressure (due to ICM B-fields and cosmic ray s 



e.g., iSkillman et al.l 120081 )) are not required, as also found by iRoncarelli et al.l (120061 ). 
Third, the nonspherical and dynamical gas/halo accretion environment present in the outer 
reaches of the numerical clusters is important in reproducing the observed properties, as is 
discussed in the next section. 



4. Discussion 

4.1. Are Galaxy Clusters in Hydrostatic Equilibrium? 

What does the particular shape of the observed and numerical cluster temperature pro- 
files imply about the dynamical state of the ICM? In particular, is this universal temperature 
profile consistent with the simplest form of hydrostatic equilibrium given in Equation (3), 
where the ICM pressure is strictly thermal? 

In Figure HI we compare the temperature profiles from the observed clusters and the 
simulations with that expected if the numerical clusters are in hydrostatic equilibrium. To 
calculate the hydrostatic temperature, we begin by making spherically-averaged 3D radial 
profiles of the dark matter density, weighted by the X-ray emission. From these radial profiles, 
we calculate the total enclosed mass as a function of radius. We then assume hydrostatic 
equilibrium to calculate the pressure derivative at each point, 

where e = p ga s/ Pdm is the ratio of gas to dark matter density and is assumed to be constant. 
We then integrate the pressure inwards, assuming P = at our outermost radial point. 
Finally, we calculate the temperature from 

firripP 

Thydro = t • (10) 

KtPdm 
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Because we normalize the temperature by the X-ray weighted average temperature between 
0.2r 2 oo < r < r 20 o, the absolute scaling of the pressure (and the value of e) does not affect 
the result. From the set of each individual profile, we then calculate the mean and standard 
deviation shown in Figure HI 

Using a x 2 -test to compare the temperatures and their standard deviations for the 
simulated clusters in Figure HJ we find that the numerical 3D radial temperature profile and 
Th y dro( r ) have a 1.6% probability of being drawn from the same parent distribution. Although 
the scatter is large, the average hydrostatic temperature profile is ~ 20% higher than the 
3D numerical clusters or observed Suzaku clusters temperature profiles. This raises concerns 
about accurately calculating cluster masses using this method, especially for applications to 
precise cosmological parameter estimation. 

In Figure 0, we investigate this issue further by plotting (M hydro (< r) — M true (< 
r))/M true (< r), where Mh y dro{< r) is the mass within a sphere with radius r calculated from 
hydrostatic equilibrium using Equation (3) and M trU e(< r) is the dark matter mass. In this 
case, we use a mass weighting to probe the physical cause of the differences in mass. Previous 
numerical simulations have shown that there is a systematic bias (underest imate) of 5-15% in 



Rasia et al. 


2006; 


Lau et al. 


2009) 



^2500 - ^500 ) (IBurns et all 12008 



this bias increases towards the edge of clusters. Beyond ~ 0.8r 2 nn, the av e rage integrated 
mass is biased low by ~15%. Recent observational work by iMahdavi et al.l (120081 ) supports 



this result in that they find hydrostatic cluster mass estimates are systematically low by 
10-20% compared to masses calculated via a weak lensing analysis out to r 50 o- Importantly, 
Figure |5] also illustrates the large scatter in the hydrostatic equilibrium calculation due to 
different dynamical states of individual clusters (even after we edited the sample for obvi- 
ous mergers). At 0.8r 20 o, cluster hydrostatic masses range from 5% overestimates to 30% 
underestimates. 
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FVR 200 

Fig. 4. — Comparison of numerical and observed cluster temperature profiles with that 
expected for hydrostatic equilibrium. The blue line and blue shading are the temperature 
distribution and standard deviations, respectively, expected if the gas in the dark matter 
potential wells of simulated halos is in hydrostatic equilibrium. The black line and grey 
shading are the X-ray weighted average temperature and standard deviations, respectively, 
for the numerical clusters. Suzaku data points are also shown for PKS 0745-191 and for two 
analysis regions (north and south) in A1795. 
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Fig. 5. — Radial profile of the fractional difference between cluster mass assuming hydrostatic 
equilibrium and the actual mass for the sample of numerical clusters. The shaded region is 
the standard deviation around the mean values. Note the increasing underestimate of mass 
with radial distance beyond ~0.3r 2 oo- 



- 16 - 



4.2. Gas Dynamics in the Outer Regions of Clusters 

Why does the intracluster medium of galaxy clusters not follow the simple prescription 
for hydrostatic equilibrium shown by Equation (3)? There are several obvious reasons which 
relate to the connection of clusters with the cosmic web and to their on-going evolution via 
accretion of gas and subclusters from filaments. 

First, unlike the view of several decades ago, galaxy clusters are not simple spheres of gas 
and dark matter that are disconnected from their surroundings. Rather, they are closely tied 
to linear filaments within the large scale structure of the Universe. Gas, galaxies, and dark 
matter are funneled along these filaments into clusters, which typically lie at the intersections 
of the filaments. As a result, accretion onto clusters is complex and nonspherical. Thus, 
azimuthal variations in VP using Equation (2) need to be folded into the calculation of 
M(< r). 

Second, this accretion also leads to bulk and turbulent gas motions in the ICM. The 
bulk gas velocities are often hundreds to even thousands of km/sec which may be detectable 
with upcoming high resolution X-ray calorimeter spectrometers (e.g., Astro-H XCS). One 
symptom of these bulk velocities and turbulent motions is the complex asymmetric temper- 
ature structures which are visible in simulated clusters as shown in Figures [T] and O The 
particular visualization in Figure was chosen to display thin isocontours of temperature, 
highlighting the various phases of the intracluster gas. This temperature image illustrates 
the complex interplay between gas of different temperatures (and velocities) accreting from 
filaments at multiple angles and thermalizing via a web of shocks throughout the cluster. 

Third, the kinetic gas motions are a significant fraction of the total energy density of the 
ICM, especially in the outer reaches of clusters. In Figure [TJ we plot the radial profiles of the 
ratios of the thermal (nkT), radial kinetic (pv%), and turbulent kinetic energy (= difference 
between total and radial kinetic energy) densities to that of the total energy density of the gas 
in our numerical clusters. To calculate the kinetic energy density of the gas, we first find the 
center of mass velocity of the halo. We then compute the difference between a cell's velocity 
and this center of mass velocity. This then gives the kinetic motion relative to the halo. Both 
the kinetic energy and thermal energy are calculated using an X-ray emission weighting. Near 
the cluster centers (< 0.2r 2 oo); the radial kinetic and turbulent kinetic energies contribute 
only « 15% to the total energy of the ICM, consistent with mild biasing of the hydrostatic 
mass estimates in the central regions of clusters. However, at r 2 oo, the radial kinetic energy 
equals the thermal energy. Furthermore, the turbulent energy density at the cluster edge 
is now ~ 50% of the thermal energy and, thus, provides significant pressure support. So, 
it is not surprising that hydrostatic equilibrium increasingly fails as a valid assumption as 
we get closer to the cluster periphery. Turbulent gas motion is an important ingredient in 
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the pressure support of galaxy clusters. For Equation (2) to result in an accurat e estimate 



of clu ster mass, it must contain a dynamical gas pressure term (~/w 2 ) (see also lLau et al. 
2009h . 
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Fig. 6. — Isocontours of temperature for a numerical cluster. The field of view is 13.85 
h~ x Mpc. The colors correspond to the following temperatures: cyan = 10 5 K, magenta = 
3 x 10 5 K, blue = 10 6 K, orange = 3 x 10 6 K, red = 1.6 x 10 7 K, and white = 5 x 10 7 K. 
This new visualization was created using a ray casting module recently developed within yt. 
This image was made by first defining a red-green-blue-alpha transfer function based on any 
grid quantity. These emissivities are then integrated along a given ray from back to front. 
Note the complex interplay in temperature between the filaments, the cluster periphery, and 
cluster core. 
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Fig. 7. — Energy density ratio profiles for the sample of numerical clusters. The thermal, 
radial kinetic, and turbulent kinetic energy densities are normalized to the total energy 
density within each shell. 
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5. Summary and Conclusions 

In comparing new Suzaku X-ray observations out to ~ r2oo f° r two rich galaxy clusters, 
PKS 0745-191 and A1795, with numerical clusters extracted from a large cosmological vol- 
ume, we find that the observed clusters follow universal temperature, density, and entropy 
profiles. This excellent agreement between observations and simulations suggests that the 
simple adiabatic gas physics of the simulations appears to adequately define the thermody- 
namics outside the cores of these clusters. With the present Suzaku X-ray data, cooling, 
non-gravitational heating (via e.g., AGNs), magnetic fields, or cosmic ray pressure effects 
are not required. 

These profiles, however, are not consistent with simple hydrostatic equilibrium in the 
outer cluster regions, suggesting that the ICM pressure support does not arise solely from the 
thermal energy density. Using cosmological numerical simulations, we show that nonspheri- 
cal accretion from large-scale filaments and the resulting gas motions (radial and turbulent 
kinetic) have a comparable energy density to the thermal energy density near r 2 oo- Assum- 
ing that thermal gas motions provides all the pressure support results in an increasingly 
inaccurate estimate of the total cluster mass with radial distance from cluster centers. 
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